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SPECTRAL INDICATIONS OF DENSITY VARIABILITY IN THE CORONA OF AD LEONIS 

A. Maggio 1 , J.-U. Ness 2 

Draft version February 2, 2008 

ABSTRACT 

Direct comparison of high-resolution X-ray spectra of the active dMe star AD Leo, observed three times 
with Chandra, shows variability of key density-sensitive lines, possibly due to flaring activity. In particular, 
a significant long-duration enhancement of the coronal density is indicated by the Ne IX and Si XIII He-like 
triplets, and possibly also by density-sensitive Fe XXI line ratios. 

Subject headings: stars: individual (AD Leo) — stars: activity — stars: coronae — stars: late-type — X-rays: 
stars 



1. INTRODUCTION 

The plasma density is a crucial parameter for understand- 
ing the structure of stellar coronae. Rather than being uni- 
formly distributed on large spatial scales, as occurs in the 
photosphere, the coronal plasma is trapped by surface mag- 
netic fields within scale lengths typically much smaller than 
the stellar radius, and is likely to be very inhomogeneous in 
density. The patchy appearance of the solar corona in X-rays 
is mainly a density effect; since the radiative emission of an 
optically thin plasma excited by collisional excitation is pro- 
portional to the square density, the observed irradiance traces 
preferentially high-density "compact" regions, confined by 
small-scale magnetic fields. The plasma density is also ex- 
pected to vary in time owing to the dynamics of these fields 
and of the plasma itself; in particular, the coronal density 
is expected to increase locally as a consequence of impul- 
sive heating events (flares) which trigger an upward motion 
of denser plasma from deeper atmospheric layers ("chromo- 
spheric evaporation"). 

The measurement of plasma densities in stellar coronae is 
a relatively recent practice in X-ray spectroscopy; Chandra 
and XMM-Newton high-resolution gratings provide us with 
a number of density-sensitive emission lines, which are now 
routinely employed to this purpose. However, density esti- 
mates and their interpretation are more difficult and challeng- 
ing than expected, for several reasons; first, even in relatively 
strong X-ray sources the relevant emission lines may be dif- 
ficult to measure owing to low ion abundances, severe blends 
with other spectral features, or insufficient instrumental sen- 
sitivity in the relevant wavelength ranges; secondly, different 
spectral diagnostics probe coronal regions at different tem- 
peratures, and they often yield density estimates which are 
difficult to reconcile; finally, the lack of spatial resolution in 
observations of stellar coronae allows us to obtain only "ef- 
fective" densities, averaged over all the coronal structures in 
the visible hemisphere and weighted by the plasma irradiance, 
so that any interpretation becomes model dependent. 

For the above reasons, it is currently easier to look for vari- 
ations with time in the plasma density in a given coronal 
source, rather than comparing absolute values derived from 
snapshots of different stars. In this respect, energetic events 
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like stellar flares are natural experiments to examine, and 
they have been the subject of several investigations in recent 
years; nonetheless, convincing evidence of density variations 
are scanty in the literature, with notable exceptions such as 
the flare recen tly observed in the nearby dMe star Proxima 
Centauri (Gud el et all2002l) . 

In this paper we present the first evidence of significant 
long-term density variations in the corona of AD Leo, an- 
other well-known dMe flare star observed several times with 
the Chandra X-ray Observatory ( Weisskopf et al. 2004). The 
strength of the present evidence comes from the direct com- 
parison of high-resolution X-ray spectra taken at different 
times, with little need of data adjustments or problematic 
atomic models. On the other hand, the interpretation of the 
results is unclear, but gives some new hints on issues under 
debate in stellar X-ray spectroscopy and coronal physics. 

2. TARGET AND OBSERVATIONS 

AD Leo is one of the X-ray brightest single dMe stars cur- 
rently known, having quite a stable quiescent X-ray luminos- 
ity, L x = 3-5 x 10 28 erg s " 1 (0.5-4.5 keV ban d) over a 17 yr 
period, as determined by Fav ata etafl fcOOOal) . On the other 
hand, it is also a quite variable X-ray source on short time 
scales, owing to its characteristic coronal magnetic activity 
(see e.g. lAudard et alJ2000l) . 

AD Leo was observed twice with the Chandra Low-Energy 
Transmission Grating Spectrometer (LETGS) in January and 
October 2000, as part of the GTO program, and once again in 
June 2002 with the High-Energy Transmission Grating Spec- 
trometers (HETGS). A detail ed analysis of the Oc tober 2000 
observation was presented by Maado et al. (2004), including 
a reconstruction of the plasma emission measure distribution 
(EMD) with temperature, and evaluation of the plasma den- 
sities at various temperatures by means of the He-like triplets 
and Fe XXI lines. Here we will focus on the comparison of 
this observation with the previous one, since they have about 
the same exposure time (48 ks and 46 ks, respectively), and 
we will also briefly discuss results based on the shorter Jan 
2000 observation (lOks). 

The data were retrieved from the Chandra archive and re- 
processed with the Chandra Interactive Analysis of Observa- 
tions (CIAO V3.1) software. We used the tool FULLGARF to 
calculate effective detector areas. From each observation we 
have extracted either the LETGS spectrum (6-140 A), or both 
the first-order Medium-Energy Grating (MEG, 1-25 A) and 
High-Energy Grating (HEG, 1-18 A) spectra, simultaneously 
provided by the HETGS; photons from the two dispersion di- 
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rections were coadded to maximize the S/N ratio. We also 
extracted the spectra in consecutive time bins in order to cre- 
ate the light curves, as described in the following section. 



TABLE 1 
Density-sensitive Fe XXI lines 
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FIG. 1. — AD Leo X-ray light curves (net observed flux in the 7-25 A wave- 
length range) derived from the Chandra observations reported in the legend. 
The fluxes were computed by integrating the dispersed spectra extracted in 
equal time bins, within each observation. 



3. ANALYSIS AND RESULTS 
3.1. X-ray light curve 

Figure [2 shows the X-ray light curves obtained from the 
three Chandra observations: dispersed spectra were extracted 
from successive time intervals of equal duration, and total X- 
ray fluxes were computed in the common wavelength range 
7-25 A, after background subtraction and correction for the 
instrument effective area 

No strong isolated flare is visible during these observations, 
but the Jan and Oct 2000 light curves both show an initial 
emission level higher than the quiescent one (/ x rj 10~ n erg 
cm" 2 s" 1 ) by a factor < 2, followed by a clear decline; this 
behavior resembles the decay phase of a flare, but we note 
that in the Oct 2000 light curve the decline is much slower 
than the expected exponentia l rate. The latter observation 
was recently studied by van den Bessel aar et al . (2003), who 
performed a separate analysis of the first 12ks segment and 
concluded that the high-temperature tail of EMD was en- 
hanced wit h respect to the fo llowing time segment; on the 
other hand, Maggio et al. (2004) argue that a rotational mod- 
ulation effect cannot be completely excluded because the du- 
ration of the October 2000 observation is a sizable fraction 
(~ 20%) of the phot ometric rotational period (f 10t = 2.7 d, 
Snie sman & Hawle\ll986l) . 

3.2. Analysis and Comparison of Spectra 

We have focused our attention on the spectral ranges in- 
cluding the Si XIII, Mg XI, Ne IX, and O VII He-like triplets, 
and on selected density-sensitive Fe XXI lines. Note that the 
first three triplets are visible in all LETGS and HETGS spec- 
tra, while the O VII triplet can be observed only with the 
LETGS and the HETGS/MEG. With regard to the Fe XXI 
lines, we have considered both 2s 2 2p3d transitions to the 
2s 2 2p 2 ground state which fall at wavelengths in the range 
12.28-12.53 A, and transitions from 2s2p 3 states which are 
visible at w avelengths 102.2-145.7 A only in LETGS spectra 
jNess et alJ2004l for more details). 
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NOTE. — Col 1: theoretical wavelength. Col 2: critical density, and 
arrow indicating the emissivity trend for increasing N e . Col 3: observed 
counts or upper limits. Col 4: LETGS effective area. Col 5: ISM ab- 
sorption factor, assuming N H = 3 X 10 18 cnT 2 <Cullv et all 19971) . Col 6: 
photon flux ratio. Col 7: estimated electron density. 



For the comparison of LETGS and HETGS spectra (Fig.|2) 
we proceeded as follows. First, we have fitted the relevant 
lines in the MEG and HEG spectra using the program CORA 
iNess & Wichmannl2 002): the smooth representation of each 
spectral region obtained with the above line fitting was con- 
verted to energy fluxes by use of the effective areas and ex- 
posure times, and finally rebinned to the resolution of the 
LETGS spectra; this procedure allows the proper overlay of 
the HETGS and LETGS data for ease of direct comparison. 

Note that the line fitting is performed with a maximum like- 
lihood method applied to the total source+background spec- 
trum. In each spectral region the strongest lines are fitted 
with normalized Moffat profiles (Lorentzians with an expo- 
nent (3 = 2.4), representing the instrumental line spread func- 
tion. The instrumental background is not subtracted from the 
total measured spectrum, but is instead added to the model 
spectrum consisting of the sum of all line emission compo- 
nents. The line counts are obtained as normalization factors 
of the lineprofiles. 

Figure |2] shows the HETGS vs. LETGS (October 2000) 
comparison for the Si XIII and Ne IX triplets; the original 
HETGS/MEG spectrum in the region of Ne IX is also shown. 
What is immediately evident is the excess emission at the po- 
sitions of the Si XIII and Ne IX intercombination lines (6.69 A 
and 13.55 A, respectively) in the LETGS spectrum with re- 
spect to the HETGS case. No such enhancement is visible in 
the O VII (21.6, 21.8, 22.1 A) and Mg XI (9.17, 9.23, 9.31 A) 
triplets (not shown); note that the latter is affected by low sig- 
nal to noise ratio because of the relatively low Mg abundance 
in the corona of AD Leo (Maggio et al. 2004). 

We have verified that the above enhancements are still vis- 
ible, although somewhat reduced, also in the last 32 ks of the 
October 2000 observations, i.e. where the X-ray emission is 
again at its quiescent level. Some enhancement is also vis- 
ible in the Ne IX intercombination line of the January 2000 
LETGS spectrum, although the shorter exposure time makes 
the signal to noise ratio too low for any robust indication. 

We have also carefully re-examined the long-wavelength 
(A > 100 A) region of the LETGS spectrum, for evidence of 
density-sensitive Fe XXI lines; in particular, we have searched 
the lines listed in Tabled and measured total counts or upper 
limits, as appropriate. We confirm the detection of the strong 
128.7 A line (Fig.[2j>, and also of weaker spectral features at A 
102.2, 121.2, and 145.65 A, already tentatively identified by 
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FIG. 2. — Upper panels: LETG Si XIII and Ne IX triplets compared with scaled HETGS line fitting models. Lower left: LETG spectral region including 
Fe XXI lines at 121.2A and 128.7 A. Lower right: Ne IX triplet in the MEG spectrum of AD Leo, with 18 best-fit line components. 



Magg io etall (|T004). Unfortunately, all the latter three fea- 
tures are characterized by quite a low S/N ratio (2-3cr), and 
there are residual uncertainties in the wavelength calibration 
of the LETGS that suggest caution in their identification; they 
fall at the nominal position of density-sensitive Fe XXI lines, 
and we will treat them as such in the following. The densities 
(or upper limits) implied by these diagnostics were estimated 
from the predicted line ratios (Fig. 0, according to APED 
(Plasma Emission Database VI. 3, Smi th et alJl2Q01l) : all the 
observed ratios, properly corrected for the instrument effec- 
tive area and ISM absorption factor, consistently indicate a 
density A^ e in the range 10 12 °-10 12,2 cm" 3 , within statistical 
uncertainties, with the exception of the 117.51/128.73 ratio. 
Note that the non-detection of the 1 17.51 A line is puzzling, 
because this line is predicted to be relatively strong and very 
weakly dependent on the electron density, hence we suspect 
that its theoretical emissivity is overestimated. 

It is more difficult to measure the Fe XXI line fluxes at short 
wavelengths because the spectral region 12.2-12.6 A is espe- 
cially crowded with hundreds of relatively faint iron lines, 
mostly from Fe XIX- XXII. In particular, we have searched 
for the density-sensitive lines at 12.284 A, and 12.327 A. 
The former has the highest emissivity at densities below 
2V e w 10 12 7 cm -3 and it is clearly present in all X-ray spec- 
tra we have analyzed, while the second line is expected to 
become enhanced only in the high-density regime, above 
N e w 10 124 cm~ 3 , and it is not convincingly detected in any 
of the spectra. 
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FIG. 3. — Fe XXI predicted density-sensitive Fe XXI line ratios, for 
the temperature range 

10 6.8_ 10 7.2 K (lower and 

upper borders). The shaded 
areas indicate the density ranges consistent with the data; the vertical lines 
are placed at the high-density limits. 



4. DISCUSSION AND CONCLUSIONS 

The measurement of coronal plasma densities by means of 
X-ray spectral diagnostics is a difficult task (see §1). For ex- 
ample, the Ne IX triplet is heavily blended with Fe XIX lines, 
and to a lesser extent with Fe XX- XXI lines. A detailed anal- 
ysis of this spectral regi on in Chan dra and XMM-Newton 
spectra of Capella (Ness et al J 120031) showed that APED is 
sufficiently accurate and complete that all significant observed 
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lines in this region can be reasonably identified in the Chan- 
c/ra/HEG spectra, but spectra at lower resolution may be in- 
adequate to derive reliable results. In fact, independent anal- 
yses of the Ne IX triplet in the sam e (October 2000) LETGS 
spectrum of AD Leo , performed bv lNess et al.l (|2002) and by 
Maggi oet alJ fc004l) . resulted in quite different estimates of 
the plasma density owing to the different treatment of blends 
by these authors. 

An appar ent inconsist e ncy em erges also by comparing the 
results in [ Maggio et al. (2004) and in the recent work by 
iNess et al.l fe004l) . who included the Chandra/HETGS ob- 
servation of AD Leo in a large sample of grating spectra 
of stellar coronae. While the former paper reports densi- 
ties N e ~ 4 x 10 11 cm" 3 at T ~ 10 6i> K (peak emissivity tem- 
perature of the Ne IX) and N e ~ 2 x 10 12 cm" 3 at T ~ 10 7 K 
(Fe XXI), no evidence of especially high densities appears in 
the survey by Ness et al. ( 20041). 

The present investigation resolves this issue; the direct 
comparison of the LETGS and HETGS spectra of AD Leo 
suggests that the star was observed in different coronal emis- 
sion regimes, a high-density state in October 2000 and a low- 
density state in June 2002. The high density is strongly sup- 
ported by the enhanced flux observed at the positions of the in- 
tercombination lines of the Ne IX and Si XIII He-like triplets; 
our investigation of density-sensitive Fe XXI lines in the Oc- 
tober 2000 LETGS spectrum appears to be consistent with the 
above conclusion, although there remain uncertainties in the 
identification of some of these lines. Note that densities of a 
few x 10 12 cm" 3 in the corona of AD Leo have already been 
reported bv ISanz Forcada & Micelal J2002I) . based on EUVE 
spectra. 

The above picture is not entirely clear for two reasons; 
first, a density increase should affect both the intercombina- 
tion (0 line and the forbidden (/) line in each He-like triplet, 
while we see only an enhancement of the i line; secondly, the 
temperature-sensitive ratios G(T) = (i+ f)/r, with r denoting 
the resonance line, is significantly higher in the October 2000 
LETGS spectrum than in the HETGS spectrum, with the con- 
sequence that, apparently, the plasma temperature was lower 
in the high-density state, a result which is difficult to under- 
stand. Since the Ne IX / line is strongly blended with Fe XIX, 
which forms at a higher temperature (T ~ 10 6 9 K), we have 
checked whether the enhancement of the i line could be ex- 
plained by a contamination effect; none of the strongest lines 
from other comparably hot ions appears to be significantly 
different in the two spectra, and hence we have no reason to 
suspect an increased Fe XIX line emission in the Oct 2000 
corona of AD Leo. Another possibility is that of a cooling 
plasma in which the population of the ls2/;> 3 P levels (and 
hence the i line) is preferentially enhanced with respect to the 
3 S level (/ line) in case of strong recombination effects; how- 
ever, the expected recombination time scales are much shorter 



than the observation length, hence a continuous replenishment 
of plasma is required to explain our finding. We conclude that 
the temperature indicated by the Ne IX G(T) ratio is not reli- 
able, for reasons yet to be understood. 

Another puzzling result comes from the comparison of the 
O VII triplet in the different observations: contrary to naive 
expectation, the density-sensitive f/i ratio turns out to be sig- 
nificantly higher in October 2000 than in June 2002 (3.2 ±0.6 
vs. 1.9 ±0.6, respectively), indicating a lower density in the 
first case (logA^ = 9.9 ±0.7 and 10.5 ± 0.3 cm" 3 , respec- 
tively). Instead, the temperature-sensitive G(T) ratios are in 
perfect agreement, and they always indicate an "effective" 
formation temperature of rj 2 x 10 6 K, i.e. coincident with the 
O VII emissivity maximum. Again, this is not easy to explain, 
because the plasma emission measure distribution vs. tem- 
perature in the corona of AD Leo peaks at rj 8 MK, and hence 
one expects an important contribution to the O VII line emis- 
sion coming from plasma at temperatures higher than 2 MK. 
This anomalous behavior of the G( T) ratios was also noted in 
the survey of coronal densities by Testa et all (120041) . 

In conclusion, our detailed comparison of the different 
Chandra observations of AD Leo indicates that the plasma 
density in the corona of this active star may change substan- 
tially. The coincidence of the observed density enhancement 
with a higher X-ray emission level is suggestive of a dynamic 
phenomenon, possibly related to flaring activity. However, the 
observed variability of the X-ray emission can be due also to 
active regions in the stellar corona coming in and out of sight, 
due to rotational modulation, or due to the natural lifetime of 
the magnetic fields which provide confinement and heating 
of the plasma. For this reason, a significant variation of the 
temperature may not occur. 

In any case, the interpretation of widely used spectral diag- 
nostics, such as the He-like triplets or Fe XXI density-sensitive 
lines is certainly complicated by the likely existence of den- 
sity inhomogeneities in the coronal plasma, besides obvious 
uncertainties in the atomic physics and in our ability to re- 
solve all the relevant spectral compone nts wi t h the available 
data. In fact, as already pointed out by Gudej ( 120041) . the ob- 
served line ratios may not describe "local densities" but rather 
the steepness of the density distribution with temperature. In 
this respect, predictions of the coronal X-ray spectra, taking 
into account physically realistic distributions of the magnetic 
structures, would be the next step in achieving progress in 
stellar coronal physics and X-ray spectroscopy. 
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